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Abstract. - We have investigated the effect of Yb substitution on the Pauli limited, heavy fermion 
superconductor, CeCoIns. Yb acts as a non-magnetic divalent substituent for Ce throughout the 
entire doping range, equivalent to hole doping on the rare earth site. We found that the upper 
critical field in (Ce,Yb)Coln5 is Pauli limited, yet the reduced (H,T) phase diagram is insensitive 
to disorder, as expected in the purely orbitally limited case. We use the Pauli limiting field, the 
superconducting condensation energy and the electronic specific heat coefficient to determine the 
Wilson ratio (Rw), the ratio of the specific heat coefficient to the Pauli susceptibility in CeCoIns. 
The method is applicable to any Pauli limited superconductor in the clean limit. 
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Introduction. — Heavy fermion (HF) systems have 
been an ideal playground for investigating unconventional 
superconductivity (SC) since the discovery of SC in this 
class of materials [HE]- Much of the attention has been 
focused on the symmetry of the superconducting order 
parameter and the interplay /competition between SC and 
magnetism [3J. Such investigations laid ground for mag- 
netism as the origin of Cooper pairing, since SC in HF 
seems to occur invariably in close proximity to a T = 
magnetic instability [3]. Their large specific heat (C) 
anomaly at the SC transition and the absence of SC in the 
non-magnetic La analogs imply that Cooper pairs form 
out of heavy quasiparticles (QP). Thus, the heavy mass 
and SC appear to originate from the same mechanism. 
Yet, the relationship between SC and the Kondo lattice 
physics, at the heart of the HF problem, has only recently 
come to spotlight. In particular, a new model of super- 
conductivity for the 115 family of heavy fermions shows 
that the composite heavy quasiparticles form only when 
the system becomes superconductor [5]. 

CeCoIn 5 is an ambient pressure SC with T c — 2.3K [6 
and has the unique feature of an antiferromagnetic (AFM) 
quantum critical point located near the upper critical field 
H c2 , indicating that AFM is superseded by SC [IE]. More- 
over, the change of the SC transition from second to first 
order for transition temperatures T c < T* ~ 0.7 T c o [9] 
combined with the discovery of a second SC phase in the 



large B/T region of the phase diagram lead to the sugges- 
tion that CeCoIn5 is the first realization [TOllTl] of a Fulde- 
Ferrell-Larkin-Ovchinnikov (FFLO) state [HE]. Subse- 
quent NMR measurements have shown that only some of 
the NMR lines broadened within the second SC phase, 
consistent with local moment ordering |14j . A more re- 
cent neutron diffraction experiment has found that this 
second SC phase carries AFM order that collapses at the 
same upper critical field at which SC is destroyed [IB] . 
CeCoIn 5 is thus the first example of magnetic order being 
stabilized by superconductivity, suggesting a ground state 
differing from the one proposed by FF and LO. 

The unconventional d x 2_ y 2 gap symmetry in CeCoIns 
has been established based on angular dependence of C 
[16] and thermal conductivity (THE], as well as point 
contact spectroscopy [TH] ■ Recently, a resonance peak has 
been discovered below T c in inelastic neutron scattering 
20 , suggesting a magnetically mediated pairing in anal- 
ogy with the high-T c cuprates. CeCoIn 5 has also been 
a model system for investigating the emergence of coher- 
ence in a Kondo lattice. A phenomenological two-fluid 
model has been successfully used to describe the crossover 
from single-ion Kondo behavior to coherent heavy fermion 
ground state in Ce\- X La^ Coins [21] . 

Here we report a specific heat investigation of CeCoIns 
as a function of Yb substitution to Ce. Yb acts as a non- 
magnetic divalent substituent for Ce throughout the en- 
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tire doping range, equivalent to hole doping on the rare 
earth site. The divergence of the specific heat in the nor- 
mal state (at H = 5T) is moderately suppressed with Yb 
doping, as a result of the Ce Kondo-lattice dilution. The 
upper critical field in (Ce,Yb)CoIn 5 appears to be Pauli 
limited, as in the pure compound, yet it exhibits a scal- 
ing expected in the purely orbital limit. To the best of 
our knowledge, this is the first report that the upper crit- 
ical field is insensitive to the amount of doping in a Pauli 
limited superconductor. Furthermore, we show how the 
Wilson ratio (Rw), the ratio of the specific heat coeffi- 
cient to the Pauli susceptibility, can be determined from 
the superconducting condensation energy in a Pauli lim- 
ited superconductor. The method yields a Wilson ratio 
consistent with the expected value in pure CeCoIn 5 . 

Crystal Growth and characterization. CeCoIns 
crystallizes in the tetragonal HoCoGas structure. Sin- 
gle crystals of Ce\- x Ln x Coins (Ln =Yb,Lu) were grown 
from excess In flux [6] . The lattice parameters were deter- 
mined from Rietveld refinement of powder X-ray diffrac- 
tion patterns, using Si standard. The effective concentra- 
tions were determined with energy dispersive X-ray spec- 
troscopy (EDS) on the measured single crystals, as well 
as proton-induced X-ray emission microprobe (PIXE) on 
a mosaic of crystals from the same batches. The mag- 
netic susceptibility was measured using a vibrating sample 
SQUID magnetometer in a field of IT or higher applied 
perpendicular to [001]. The heat capacity was measured 
using the standard adiabatic heat pulse technique in a 
3 He-PPMS. The resistivity (p) was measured using the 
standard four wire technique with an AC resistance bridge. 
Low resistance contacts were obtained by spot- welding Au 
wires. 

Fig. Q] shows the doping evolution of characteristic pa- 
rameters in Cei_ x Lnj;CoIn5 (Ln =Yb,Lu). The effec- 
tive Yb concentrations, as determined with either EDS 
or PIXE, are close to the nominal values for small x nom 
but show large distribution around x nom = 0.5, as in- 
dicated by the error bars in Fig. QJi. Phase separation 
between pure YbCoIn 5 and Ce 1 _ 2 .Yb 2; CoIn5 is likely the 
reason why we could not reach effectively xyt > 0.3. In 
fact, for x nom > 0.7, the batches yield essentially pure 
YbCoIns, with very few Yb-substituted CeCoIns crystals. 
The difference between EDS and PIXE values reflects this 
difference between a single crystal and the average concen- 
tration of the mosaic of crystals. For simplicity, nominal 
concentrations will be used in the rest of the paper. The 
lattice volume decreases systematically with Yb and Lu 
doping due to the lanthanide contraction (see Fig. [lb). 
The Curie- Weiss moment (J-eff (per formula unit), ob- 
tained from linear fits to inverse magnetic susceptibility 
is suppressed below the Ce 3+ moment (2.54/ie) with Yb 
as with Lu doping (see Fig. QJi). This indicates that Yb 
substitutes as a non-magnetic Yb 2+ ion, resulting in a di- 
lution of the Ce lattice. As such, the Yb doping is not an 
isoelectronic substitution, unlike La or Lu, but is equiv- 
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Fig. 1: (Color online) (a) Effective (xyt) vs nominal (x nom ) 
concentrations of Yb in Cei- X Yb x Coins, as determined from 
EDS and PIXE. (b) Critical temperature T c vs x nom in Cei-^ 
YbxCohis, as determined from magnetic susceptibility (x), 
specific heat (C) and resistivity (p). (c) Lattice volume vs 
Xnom in Cei-z Ln x Coins (Ln =Yb, Lu), as determined from 
powder X-ray diffraction, (d) Effective Curie-Weiss moment 
fi e ff (in units of Bohr magneton) vs x nom in Cei-^ Ln^CoIns 
(Ln =Yb, Lu), as determined from \. 



alent to hole doping. The absence of Curie- Weiss behav- 
ior in pure YbCoIns and its small Sommerfeld coefficient 
(~ 11m J / K 2 mol) shows that it is not a heavy fermion. 

Results and Discussion. — Figures and [2b show 
the zero field superconducting transition in p and C for 
various Yb (EDS) concentrations. The transitions in p are 
sharp, except for x — 0.6 which shows a double jump, in- 
dicative of an inhomogeneous sample, consistent with the 
large error bars found in the EDS (see Fig.[T^,). Two differ- 
ent crystals have been measured in p from most batches 
and they exhibit very similar T c 's, determined from the 
onset of non-zero resistance, and shown in Fig. [Tb>. A 
fairly sharp, A— like anomaly is observed in C, measured 
on the same single crystals for x — 0.1,0.2 and 0.3, also 
shown in Fig. [2b>. The SC anomaly is smaller and broader 
for x = 0.55 and 0.8. The onset of the jump in ^ de- 
fines T c for all samples. Two samples have been measured 
for x = 0.55 and both exhibit a broad jump (but with 
similar T c 's) possibly due to the doping inhomogeneity. 
Overall, the resistive T c is in good agreement with the 
thermodynamic T c determined from y and x, as shown 
in Figure LTb, except for x = 0.8 with T§ > T c c . The 
ratio -=Sr has been determined with 70 taken as the % 
value linearly extrapolated to T = in the normal state 
at H = 5T. In the standard BCS theory, this ratio is 
expected to be 1.43 in the weak coupling regime. For Yb 
concentrations x — 0.1, 0.2 and 0.3, the SC appears to be 
bulk, in contrast to x — 0.55 and 0.8 with significantly 
smaller -==- values. 

70 T c 

The temperature and field dependence of the electronic 
specific heat £ in Yb-doped CeCoIn 5 is shown in Fig. [31 
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Fig. 2: (Color online) (a) p vs T in single crystals of 
Cei-zYb^CoIns. (b) Electronic specific heat | vs T in the 
same crystals at H — 0. The indicated concentrations are 
nominal. 

for XYb = 0.1,0.3,0.55 and 0.8 (nominal). The electronic 
contribution is obtained by subtracting the phonon contri- 
bution, determined from C of YbCoIns. £ has a divergent 
T— dependence down to T c2 (H) for all Yb concentrations, 
with little change in slope with increasing magnetic field. 
Moreover, the divergence of Sr extends down to lowest 
T ~ 0.5K at H = 5T in these samples, a characteristic 
shared with the pure CeCoIn 5 . This divergence has been 
attributed to a field tuned QCP near H® 2 in CeCoIn 5 [7115] , 
As seen in Fig. [2p and Fig. [3j the divergence of M becomes 
weaker with increasing Yb concentration and the corre- 
sponding 70, obtained from linear extrapolation of S vs T 
at H — 5T, decreases systematically with xyb, consistent 
with a dilution of the Kondo lattice (see table [1]). 

Figure [3] also shows the magnetic field suppression of 
the superconducting transition, for magnetic field applied 
parallel to the [001] axis, with a broader and smaller jump 
as the field is increased, due to the presence of vortices. 
Unlike in the pure case, the SC transition in the Yb-doped 
compounds does not sharpen into a first order anomaly in 
the T=0 limit, as seen in Fig [3J likely due to disorder 
introduced by doping [22] , The corresponding H-T phase 
diagrams (deduced from C) are shown in Fig. 0k and Hb. 
The upper critical field (H C 2(T)) is determined from the 
kink in the entropy S, corresponding to the midpoint of 
the specific heat jump seen in Fig. [3J S is obtained by 
integrating £, following a quadratic |23j extrapolation of 
C vs T to T = 0. For xyb = 0.55, the average critical 
field of two samples is shown. In pure CeCoIns the first 
order SC transition [9] as well as the H C 2{T) fits [24] give 
clear evidence for a Pauli limited H C 2 . Despite the absence 
of a first order transition in Yb doped compounds, the 
Pauli limit [2 5) still applies: the extrapolation based on 
the standard expression [26], H° 2 — —0-7^fr-T c , yields an 
orbital critical field far in excess of the observed transition 
field, see table Q] 

Figure [4^ shows that the H C 2 data for all Yb concentra- 
tions collapse onto a single curve when scaled by the initial 
slope at T c . Such a scaling is expected in the purely orbital 
limit, as the reduced critical field vs the reduced temper- 
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Fig. 3: (Color online) Electronic specific heat ^ vs temper- 
ature for nominal Yb concentrations of (a) xyb = 0.1, (b) 
xyb ~ 0.3, (c) Xyb = 0.55 and (d) xyt — 0.8, and with mag- 
netic fields ranging from to 5 T, applied parallel to [001]. 

ature -y- is only weakly dependent on the disorder level 
in this limit [26) . In the Pauli limit there is no reason to 
expect a similar scaling, however the orbital mechanism 
is still predominant in the zero field limit. The implica- 
tion of this scaling is that the Maki parameter [22], the 
ratio of the orbital critical field to the Pauli limiting field 

• a = ^Hp^ ' * s independent of xyb- This is not a trivial 
result, knowing that a decreases under pressure [27] . a 
was estimated to be 3.6 for £f||[001] in pure CeCoIn 5 [5]. 
We have used this value to determine Hp in the Yb-doped 
samples from H® 2 , see tabled]). 

In a d-wave BCS superconductor, the superconducting 
condensation energy, U c , is related to the specific heat 
jump AC at the SC transition via the standard relation 
[28] : U c = ( A( p Tc . Since £ in the normal state is di- 
vergent in CeCoIn 5 and since it is likely a strong coupling 
superconductor [B], the use of this formula is questionable. 
Alternatively, we have determined U c directly from inte- 
gration of ^ a {S n — S s )dT up to T c , following extrapolation 
of the £ data to T=0. We have excluded the x=0.55 and 
0.8 samples from the analysis since the lowest temperature 
(0.4K) in the data does not allow a proper extrapolation 
to T=0. The results are listed in table [T] The obtained 
value of U c for pure CeCoIns (1.43 J mol -1 ) compares 
favorably with the value (1.34 J mol -1 ) determined from 
magnetization [29]. For the doped samples, we found that 
the condensation energy decreases with increasing Yb con- 
centration, a trend similar to the one reported for La (not 
shown) or Sn doping [30] . 

The combination of U c and Hp then allows the determi- 
nation of the Pauli paramagnetic susceptibility (xo)- In 
fact, Hp is related to U c via 25 : U c = \xoHp. This 
expression is originally derived for an s-wave supercon- 
ductor [3S] and overlooks the possibility of a finite Pauli 
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Table 1: Doping dependence of characteristic parameters in 
Cei-^Yb^CoIns. xyb '■ nominal Yb concentration, T c : critical 
temperature at H — from heat capacity, 70: the T = 
linear extrapolation of = at H — 5T in J/K 2 mol, U c : the 
superconducting condensation energy in J/mol, H° 2 (T) and 
Hp(T): the orbital upper critical field and Pauli limiting field. 
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Fig. 4: (Color online) H-T phase diagram in Cei-aiYbzCoIns. 
(a) Reduced critical field -iw 2 — vs reduced critical temperature 

~dT Ta 

Tjl^- for various Yb concentrations (nominal), (b) Upper critical 
field H C 2 vs Temperature. 



susceptibility at T=0 in the superconducting state j31) . 
For a d-wave superconductor, it only remains valid in 
the clean limit. In the presence of nodal quasiparticles 
in a d-wave superconductor, there should be a contri- 
bution to the free energy in the superconducting state 
and thus the Pauli critical field should be derived from: 
F n — XnHp = F s — XsHp, where the subscripts "n" and 
"s" refer to the normal and the superconducting states. 



This is equivalent to AF = F n — F s = (1 — ^)XnHp and 
thus U c — f (1 — —)XnHp. It is straightforward to as- 
sume that the fraction of excited nodal quasiparticles, 
should be proportional to the Yb concentration, but since 
we do not have an exact determination of this ratio, we 
will not pursue the analysis in the Yb doped compounds. 
Instead, we focus on pure CeCoIns, which is in the clean 
limit, implying that ^■■Cl. For CeCoIn 5 , xo = 10~ 4 /47r 
determined from U c via this formula is close to the value 
of the c-axis susceptibility [6], x = 1-54 10- 4 /4tt at 1.8 K. 
The resulting Wilson ratio is Rw = ^ =0.76 R w , where 
R w is the free electron value. In this method, the error 
on Rw essentially comes from the uncertainty on U c (via 
Xo)- 

The Wilson ratio of a free electron gas is defined as the 
ratio of the Pauli paramagnetic susceptibility to the elec- 
tronic specific heat (Sommerfeld) coefficient and is a uni- 

— Despite the large mass 



versal number: R w = * 



renormalization, heavy fermion systems exhibit Wilson ra- 
tios very close to the free electron gas value. This is due to 
the fact that the mass renormalization corresponds to an 
enhanced density of states at the Fermi level, and the lat- 
ter determines both the paramagnetic susceptibility and 
the specific heat coefficient [32 . The difficulty in estimat- 
ing the Wilson ratio in heavy fermion systems is associated 



with the experimental determination of the low tempera- 
ture Pauli susceptibility. In fact, the magnetic suscep- 
tibility is overwhelmingly dominated by the Curie- Weiss 
contribution of Ce moments. Here we show that this dif- 
ficulty can be overcome in the case of Pauli limited heavy 
fermion superconductors, making use of the superconduct- 
ing condensation energy. The Rw we obtained from this 
method in pure CeCoIn 5 is close to but somewhat lower 
than the expected value of Rw = 2J 2 ^ 1 R-w = ^ or 
a system of J=5/2 local moments (corresponding to Ce), 
confirming the validity of the method. 



Conclusion. — In conclusion, we have investigated 
the effect of Yb substitution on the superconducting and 
heavy fermion properties of CeCoIns . Our findings can be 
summarized as follows: (i) the suppression of the conden- 
sation energy and Sommerfeld coefficient 70 with Yb dop- 
ing is due to the dilution of the dense Ce Kondo lattice, 
(ii) the upper critical field exhibits a scaling that is un- 
usual for Pauli-limited superconductors, implying a dop- 
ing independent Maki parameter. Moreover, we introduce 
a new method for the determination of the Wilson ratio 
from the superconducting condensation energy, which is 
valid for any clean superconductor in the Pauli limit. The 
value we estimate is consistent with the expected value for 
pure CeCoIns. 
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